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We report on our theoretical investigation of the effects of the confining potential profile and 
sample size on the electron velocity distribution in (narrow) quantum-Hall systems. The electrostatic 
properties of the electron system are obtained by the Thomas-Fermi-Poisson nonlinear screening 
theory. The electron velocity distribution as a function of the lateral coordinate is obtained from 
the slope of the screened potential at the Fermi level and within the incompressible strips (ISs). We 
compare our findings with the recent experiments. 
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INTRODUCTION 



In the early electro-optical measurements performed 
on the two dimensional electron systems (2DESs), the 
electrostatic potential across the Hall bar was shown to 
exhibit local dips as a function of the lateral coordinate 
across the samplei^. The positions of these local po- 
tential variations are strongly dependent on the applied 
perpendicular magnetic field. They coincide with the po- 
sitions of the strips with finite width corresponding to 
integer local filling factors where the longitudinal con- 
ductance vanishes, i.e. <Je{x) — a XK (x) — o- yy (x) = . It 
was concluded that the "expected" quasi one dimensional 
edge states can be as large as hundreds of micrometers, 
where <Ji{x) ^ 0. It was also shown that the current can 
flow from the bulk in the magnetic field interval where 
the Hall resistance does not assume its quantized value. 
In this regime, however, the outermost edge states are 
reported to be "invisible"—. On the other hand, the lo- 
cal probe of the electrostatic potential and the longitu- 
dinal resistivity measurements^!^ using scanning force 
microscopy and single-electron transistor (SET) has in- 
dicated that the current is confined within finite regions 
across the sample. These regions were later suggested to 
be the "incompressible" regions, namely regions of inte- 
ger local Landau filling factors which are distributed in an 
inhomogeneous manner over the sample due to the elec- 
tronic nonlinear screening as well as the boundary effects 
as previously predicted^. These experiments are well 
explained by the recent theoretical work s 10 ' 11 ! 12 , which 
take into account interaction effects by exploiting the 
smooth confining potential profile within the Thomas- 
Fermi approximation and also incorporating a local ver- 
sion of the Ohm's conductivity model. These models 
contribute not only to the understanding of the induced 
electric field and current distribution, but also to the high 
precision nature of the low-temperature integer quantized 



Hall (QH) plateaus in narrow Hall bars as a function of 
the continuous lateral sample coordinate. 

Recently, the edge profile a InP/InGaAs Hall sample 
was probed in the "surface photo-voltage (SPV) spec- 
troscopy" measurements and it was found that the elec- 
tron velocity at the edges increases with increasing mag- 
netic field (B) respecting a square root behavior—, i.e. 
v c \(B) (x B 1 / 2 . In the interpretation of the data these 
authors used a model in which the electrostatic bending 
of the Landau levels (resulting from the confinement po- 
tential) was not taken into account, the electron Hartree 
potential was neglected and instead, the velocity distri- 
bution was modelled by an homogeneous induced elec- 
tric field. It was stated therein that, these measure- 
ments should be re-interpreted in the light of a "screening 
model" . It is the aim of the present paper to show that 
their original interpretation cannot be corrected even by 
including a linear screening model (see Sec. IIII[) . In par- 
allel to these developments, the nonlinear screening was 
also promoted by the importance of the recent electronic- 
Mach-Zehnder interferometer experiment o 14 ' 15 where the 
role of the electron-electron interaction on the B field de- 
pendence of the edge fields was emphasized. In these lat- 
ter experiments the electron phase deduced from roughly 
assuming a constant group velocity of v g ~ 2 — 5.1 6 
cm/sec disagrees with the single particle picture, and 
the authors argued this in favor of a possible breakdown 
of the single particle picture and the Landauer-Buttikcr 
conductance formalism. 

In the present work, we systematically analyze the 
electrostatic edge profile of narrow Hall bar samples 
using a self-consistent Thomas-Fermi-Poisson approach 
(SCTFPA) under QH conditions. In Sec. HQ we intro- 
duce a model which incorporates a constant donor den- 
sity ala the Chklovskii geometry^, and secondly (Sec lIII[) 
by a nonuniform distribution of donors. We then find, in 
Sec. IIIII and Sec. HVI the electron velocity v y , considering 
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different models, across the sample in the current direc- 
tion by 



dk v 



(1) 



where Ex.n is the eigenenergy of the single particle 
Hamiltonian with X — hk y /eB, denoting the center co- 
ordinate, and n, the Landau Level (LL) index. Here, k y 
is the conserved electron momentum in y-direction, e is 
the electron charge and B represents the strength of the 
perpendicular magnetic field. We then investigate the 
dependence of the electron velocities on the B field con- 
sidering two edge state models in Sec. IIV1 The widths of 
the ISs depending on the sample properties are examined 
in Sec. [V] We observe that the electron transport is con- 
fined within the ISs where the electron velocity decreases 
with increasing magnetic field as B^ 1 / 2 . On the other 
hand, if the center filling factor, f (0), is smaller than its 
minimum integer value 2 (since, we do not resolve the 
spin degeneracy), all the current is spread over the sam- 
ple, suggesting that the slope of the screened potential 
should be calculated at the Fermi level. We close our 
discussion with a summary section. 



II. 



THE MODEL 



The 2DES, described by the electron number density 
n e i(x), is considered to be in the x — y plane with a lateral 
confinement \x\ < b at z = and assuming translation 
invariance in the y-direction. The ionized donors also re- 
side in this plane, with the average number density no 
confined into the interval \x\ < d, where d is the sample 
width and (d — b) the depletion length with b < d. Elec- 
trostatic self consistent solution is then independent of 
the y coordinate and from the solution of the Poisson's 
equation with the appropriate boundary conditions, i.e. 
V(x — ±d, z — 0) = 0, we obtain the Hartree potential 
energy of an electron in the plane of the 2DES as 



2e 2 f d 

Vn(x) = —r- I dx' K (x, x') n e i(x'), 
« J-d 



(2) 



with k being the dielectric constant of the material and 
the electrostatic kernel^ is 



K(x,x') = In 



^(d 2 - x 2 )(d 2 ~ x' 2 ) +d 2 - x'x 



(x — x')d 



(3) 



The potential energy of an electron in the confinement 
region generated by the donors reads: 



V hg (x) = -Eoy/1 - (x/d) 2 , E = 2ne 2 nod/R, : (4) 

which can be found from Eq. {2]) using the kernel given in 
Eq. (PJJ) and replacing n e \{x') by — n . We write the total 
potential energy of an electron as V(x) = Vr(x) + Vb g (x). 



The electron number density is calculated numerically, 
within the Thomas-Fermi Approximation (TEA) 

n e i(x) = J dED(E)f([E + V(x)-fi\/k B T) (5) 



with D(E) describing the (collision-broadened) Landau 
density of states (DOS), /(a) = (1 + e a )~\ the Fermi 
distribution function and /x the electrochemical potential. 
Here, ks and T represent the Boltzmann constant and 
the electron temperature respectively. We also assume 
the electron spin degeneracy. Eq.'s {2} and (JSJ) complete 
the self-consistent schem o 10 ' 11 , which can be solved by a 
numerical iteration. For accurate convergence, we first 
perform calculations at T = and B = 0, then increase 
T at an elevated B strength, and reduce the tempera- 
ture stepwise until the relevant temperature is achieved. 
In the next sections, we first consider two distribution 
functions for the donor number density, and investigate 
the electron velocity dependence on the magnetic field, 
temperature and confining potential profile. 



III. DONOR DISTRIBUTION 

During the last decades, several boundary conditions 
were considered, ranging from infinite hard-wall poten- 
tials^ to smooth potentials ^ 16 ' 18 in order to theoreti- 
cally investigate the 2DES under QH conditions. For 
relatively large samples (d > 1 mm) the edge effects were 
considered to be dominated by the localization, thus the 
positions of the ISs were mainly predicted by the disor- 
der potential. On the other hand, for narrow samples 
(d < 15/im) the ISs were considered to be formed due to 
the electrostatic boundary conditions at the edges. Re- 
cently it was experimentally shown that, the steep po- 
tential at the edge of the sample prohibits the formation 
of the ISs and the Chklovskii picture is no longer appli- 
cable^. These results coincide with an early theoreti- 
cal calculation, based on Hartree approximation, given 
in Ref. where the edge potential is taken to be 

an infinite wall, for which ISs were not observed. On 
the other hand, if the edge profile is smooth, several in- 
compressible regions can be observed theoretically — i 18 i 20 
which are confirmed experimentally 21 . In the intermedi- 
ate case, corresponding to narrow samples, only a single 
incompressible edge strip was reported^ which was then 
supported by subsequent theoretical work o 11 ' 22 . 

In this section we consider narrow samples (lpm < 
d < 5pm) and vary the donor distribution, to investi- 
gate the widths of the ISs depending on the magnetic 
field strength. The selected donor profiles can be real- 
ized experimentally either by the uncontrollable etching 
processes or by gradually doping the sample. 

In Fig. [1] we show the two selected donor distributions 
(upper panel) and the corresponding confinement poten- 
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tials (lower panel) generated by 
Pi(u) = 



l-(u+c) 2 + (c-l) 2 ]n c 
l^W 2 



-(u~c) 2 + (l-c) 2 ]n c 
(c-l)= 



, — 1 < u < —c: 

—c < u < c; 
, c < u < 1. 



(6) 



and 



(l-c) ' S u < -c, 
P2(") = ^ «c, -c < u < c; 

c<«<1- 



(7) 



where it = x/d, and n c is a constant density preserving 
the total number of the donors in the sample. The steep- 
ness of the confinement is controlled with the dimension- 
less parameter c. In the figures, the potential energies are 
also normalized with the pinch-off energy (Eg) of the con- 
stant donor distribution. To make a connection between 
the experimental realization of such donor distributions, 
we point that, during the chemical etching in the z— di- 
rection the reaction also takes place in the x — y plane. 
Hence, the donor layer is not necessarily etched com- 
pletely at the edges and a distribution similar to p\ (u) is 
expected. Meanwhile, during the growth process of the 
wafer, donors can be distributed (in a controlled way) 
similar to, P2{u). In our calculations in both cases, we 
keep the average donor number density constant. It is 
clearly seen in Fig. [T] that the steepness profiles close to 
the edges of the sample is different for two distribution 
functions, whereas the minima of the confining potentials 
change linearly with c. As a result, different behaviors 
can be identified for the screened potentials (even with- 
out exploiting the magnetic field) arising from the mo- 
mentum (q) dependence of the Thomas-Fermi dielectric 



function, i.e. e(q) = 1 



dominated by q = 2~k / a. 



The relation between the screened potential, V scr (q) 
and the external potential V cx t(l) is given by, 

V SCT (q) = V CKt (q)/e(q). (8) 

Different steepness values lead to different dominating q 
regions which then render different characteristic screen- 
ing properties. From the inset of Fig. [T] we conclude that, 
the steepness of the potential increases much faster for pi , 
better simulating the edge profile than the doped profile. 
This difference becomes crucial, when a magnetic field is 
applied and the ISs are formed at the edges of the sample. 

In the next step we consider the effect of the elec- 
tronic screening. We calculate the screened potential 
self-consistently by solving Eq.'s J2]) and (JSJ) at T = 
and B — 0. In this limit, Eq. (J5]) is reduced to 



n el (x) = D (E° 



V(x))Q{E° 



V(x)), 



(9) 



which then becomes a linear relation between the po- 
tential and the electron distribution within the linear 
screening regime. In Fig. [5] the calculated self-consistent 




FIG. 1: The cross-section of the donor layer considering (a) 
pi (u) and (b) p2 (u) for various values of the steepness param- 
eter (0.5 < c < 1), together with the calculated background 
potential profiles (lower panels, (c) and (d)). The thin solid 
line represents a constant donor distribution (c = 1), whereas 
thick solid line corresponds to c = 0.5. The line code denotes 
a gradual increase of c with a step of ten percent. The sample 
width d and the depletion length, b are fixed and set to be 
3/Ltm, 300nm, respectively. In both cases, the donor number 
density is kept constant and chosen to be 4 • 10 11 cm -2 . The 
inset depicts the variation of the background potential calcu- 
lated at the center of the sample for pi (thin solid line) and 
p2 (broken line). 



yj 




FIG. 2: The screened potentials obtained from the bare con- 
finement potentials shown in Fig.[T] at B = and T = (thick 
lines). Also the Fermi energy for vanishing field and temper- 
ature (thin horizontal lines). The inset depicts the variation 
of Ep versus the steepness considering pi (solid line) and p2 
(broken line). 
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FIG. 3: The sample width dependence of the screened po- 
tential calculated at the center. The line code depicts the 
selected values of c and two distribution functions. 
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FIG. 4: The numerical derivative of the screened potentials 
for different steepness values calculated at the edge of the 
sample for inverse parabolic donor distribution, p%. Horizon- 
tal axis essentially presents the Fermi level, i.e. increasing of 
the u corresponds to increasing of the average electron den- 
sity. 



potentials are depicted for the considered donor distri- 
butions together with the variation of the Fermi energy 
as a function of c. On one hand, the screened poten- 
tial within the sample coincides with the E F for both 
models in © and 0. On the other hand these quan- 
tities differ strongly for both distributions, due to the 
non-linear screening pronounced above. Since the con- 
stant part of the donor distribution (q — component) 
is strongly screened, the change of the potential depend- 
ing on steepness is less pronounced for p%(x), meanwhile 
the sharp transition is observed at p\(x). This implies 
that more q components being involved in the screening 
there. This behavior can be seen from the slope shown in 
the inset of Fig. [1] In other words, less q components con- 
tribute, for pi(x), to the screened potential in the bulk, 
whereas, more q components are involved close to the 
edge of the sample. Thus the minimum of the screened 
potential changes faster than that of p%(x). The Fermi 
energies, show a similar behavior, up to a factor, which 
indicates that the average number of electrons decreases 
faster for pi (x) although the density of the donors is kept 
constant. 

The effect of the sample width on the potential pro- 
file affects the variation of the screened potential. In 
Fig. [3] the self-consistent potential at the center is plot- 
ted against the (half) sample width for the two donor 
distributions in Eq.'s ([6]) and ([7]) for the selected steep- 
ness parameters. For large samples (2d > 6/im), the 
variation of the central value of the potential is not sen- 
sitive to the steepness, since the electrons at the bulk can 
screen perfectly the confinement potential at the edges. 
On the other hand, steepness is expectedly important for 
narrower samples. This observation clarifies the domi- 
nating role of the edge profile on the electron velocity for 



narrow samples and shows that the interaction effects be- 
come important in the Mach-Zehnder type experiments, 
where the dimensions of the samples used are usually less 
than a few micrometers (2d < 3 — 4/im). 

Before proceeding with the investigation of the effects 
of the magnetic field, i.e. considering the effects of the 
IS in the presence of nonlinear electronic screening, we 
briefly discuss the slope of the self-consistent Hartree po- 
tential by only taking into account the LL quantization. 
As the magnetic field is changed, the LLs and their sep- 
arations are shifted on the energy axis. A qualitative 
understanding of the B dependence of the induced elec- 
tric field (F) within the sample can then be achieved by 
analyzing this shift. Then F(B) can be compared with 
the experimental results obtained by Ref. [lU . Our main 
argument is, even without taking the ISs into account, 
one should be able to observe the predicted behavior of 
the (average) electron velocity (v e i) at the edges of the 
sample. 

We use the derivative of the screened potential with re- 
spect to u at the Fermi level to infer v c \. This derivative 
of Vscx (u) is shown in Fig. [4] for selected values of c and 
using p2(u) as the donor distribution. We observe two 
characteristic behaviors. In Fig.0^, in the electron dense 
region (|u| < 0.9 where V SCI (u)/E < Ep/E ), the deriva- 
tive grows in the positive x direction rather slowly com- 
pared to that in the electron depleted region, as shown in 
Fig-EJx Screening is strong in the electron dense interior 
where the total potential is flat. Approaching the edge of 
the populated region, the number of electrons decreases 
and screening becomes poor. In the depleted region, the 
confinement potential is screened very poorly, thus the 
variation of the total potential there is large yielding a 
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larger derivative. For a given Fermi energy, decreasing 
the magnetic field corresponds to sweeping the x axis by 
which the electron velocity along y can be deduced. In 
the inset of Fig. [5J v e \ oc B~ x l 2 behavior is clearly ob- 
served for all c values, however, the exact quantitative 
values depend on the steepness of the edge profile. We 
observe that, for c < 0.8 the change in the electron veloc- 
ity is much more rapid in this case than the shallow edge 
profiles, indicating a strong relation between the edge 
profile and the magnetic field dependence of the elec- 
tron velocity. Combining v c \ oc B^ 1 / 2 and the calculated 
slopes demonstrate that, in the mentioned experiments^ 3 - 
the confinement is relatively steep; which was concluded 
by these authors to be the opposite. In connection, here 
we would like to stress another experiment where a simi- 
lar geometry reported in Ref. [l9j was considered. In this 
work it is ruled out that, if a negatively charged gate is 
placed on the side perpendicular to the 2DES, (in the 
experimental setup this gate is an other 2DES, obtained 
by a cleaved edge overgrowth technique) creating a sharp 
potential profile at this edge, no ISs are observed. Sim- 
ilarly, a side (gate) electrode is used to detect the SPV 
signal in the experimental setup of Karmakar et.ali 3 - and 
their conclusion contradicts strongly with the findings of 
Huber et. They also contradict with the velocity 

dependence, which we discuss next in more detail now 
also including the incompressible regions. It is clear that 
we will work with those geometries where the edge pro- 
file is neither very steep (such as an infinite wall or a 
perpendicular side gate) nor very shallow so that many 
incompressible regions can be observed at a given mag- 
netic field, within Thomas-Fermi approximation. 

In the SPV work, the Hamiltonian of the system was 
given by 



H 



2m: 



■(p - eA) 2 + eFx 



(10) 



which includes a constant electric field (F) along the pos- 
itive x axis pointing to the edge. Here m* is the effective 
electron mass and p and A are the canonical electron 
momentum and the vector potential respectively. Using 
the Landau gauge, the energy dispersion is found to be 2 ^ 



E,, 



(n + l/2)hu c 



-{F/B){X/lD 
-(m c */2)(F/5) 2 , (11) 



where E g is the energy band gap and l^, = \J h/muj c the 
magnetic length. These authors concluded that, in or- 
der to obtain the B = value and also to match the 
experimental results (see Fig. 4 of Ref. [ ll3T)]. one should 
assume that F oc B 3 / 2 . First of all, one remark is that, 
in the B = limit energy dispersion given in equation 
(|lip becomes meaningless. Secondly, in the limit of high 
magnetic field, assuming F oc B 3 ! 2 or F oc B 1 ! 2 essen- 
tially leads to similar linear behavior at the measured 
B values, as shown in the inset of Fig. [4] In the SPV 
experiments no low field [B < 1 T) measurement were 
performed, therefore we conclude that their conclusion 
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FIG. 5: The slopes of the screened potential calculated at 
the Fermi level (left panel) and within the IS (right panel), 
considering characteristic (half) sample widths of d = 1/im 
(upper panel), d = 2/im (middle panel) and d — 5/im (lower 
panel). The electron temperatures are chosen to be T — IK 
(thick lines) and T = 5K (thin lines). . 



about F oc £?V 2 j n R e f j s no t unique. Moreover 
our calculations ascertain that, even in the absence of 
ISs, the electron velocity is an inverse square root func- 
tion of the magnetic field, namely v c i oc B -1 / 2 , and as a 
consequence, F oc B 1 / 2 . 

Our simple self-consistent calculations, assuming that 
the effects of ISs are negligible, agree qualitatively well 
with the experimental findings. We also point that the 
functional form of the electric field and the interpreta- 
tion of the steepness of the potential strongly differ from 
Ref. [llj]. First of all, it is experimentally-^ and thco- 
reticallySiii shown that in the presence of a side gate, 
perpendicular to the 2DES, (simulating a hard wall po- 
tential or surface charges) the potential at the edge is 
steep. Secondly, the proclaimed B dependence of the 
electric field at the edge is not unique and we claim that 
FozB 1 ' 2 . 

The discussion above should also be reconsidered in 
the presence of IS. In the next section we do that by 
examining both the potential slope at the Fermi level 
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FIG. 6: The sample width dependence of the IS thickness 
for v (x) = 2. Calculations are performed at IK (thick lines) 
and 5K (thin lines) for three characteristic steepness values 
considering etched (left panel) and functionally doped (right 
panel) samples. Widths of the samples are selected to be 
d — lfim (top) d = 3/im (middle) and d = 5^im (bottom), 
whereas the electron depleted strips are fixed to be the ten 
percent of d. 

and at the position of the ISs to obtain a more realistic 
comparison between our theory and the recent Mach- 
Zehnder interferometry experiment o 14 i 15 . 

IV. COMPARISON OF THE E-FIELD AT 
DIFFERENT EDGE STATE MODELS 

In principle, the electron velocity at the edge of the 
2DES or the electric field at the depleted region is not 
directly measured in the SPV experiments, instead the 
slope of the potential profile is investigated as a function 
of the magnetic field deduced by the energy dependence 
of the SPV signal. On the other hand, in our calculations 
we explicitly obtain the self-consistent potential and by 
taking the derivative of the energy dispersion we can di- 
rectly calculate the electron velocity. In the previous sec- 
tion, by making use of the Thomas- Fermi approximation, 
we obtained the screened potential and claimed that the 
center coordinate dependent dispersion is given by, 

E n (X)=E n + V s l= - B =°(X). (12) 

The next step is to calculate the full V^ c f°' B ^{X) and 
to investigate its slope as a function of B,T,c,d as well 
as the long-range part of the disorder potential. 

As a standard technique£ ' 22 i 24 ' 25 , we simulate the po- 
tential fluctuations generated by the disorder by impos- 
ing a modulation potential 25 of the type 

V m (x) = V cos(fc A a;),with k x = {\+l/2)x/d (13) 



as an additive contribution to the confinement potential. 
Here Vo is the modulation amplitude and A is an integer 
to preserve the boundary conditions. 

In Fig. we show the numerical derivative of V SCT (x) 
at the chemical potential. Note that at T ^ the Fermi 
energy is no longer equal to the chemical potential, and 
it has to be calculated for the given set of physical pa- 
rameters. Here we consider A = 5. The amplitude of 
the modulation is set such that, after screening, the po- 
tential variation is at the order of %5 — 25 of Ep. As 
a rough estimate we find that Vq reduces by three or- 
ders of magnitude due to the dielectric screening, for 
GaAs, k = 12.4, and electronic screening (e(q) ~ 41, 
see e.g. the expression given in the caption of Fig. 2 of 
Ref. [9( and the related text). Regardless of the variation 
in d, T and Vo, the slope of the screened potential ob- 
tained at the chemical potential exhibits the previously 
observed F cx B^ 1 / 2 form; except the case where the 
modulation is so strong that slope remains unaffected 
at d = 1/zm with Vq = 300 meV. The hint to under- 
stand this exception is found in Fig. [5b, where we show 
the slope calculated inside the IS at the position corre- 
sponding to v(x) = 2. We see that the derivative of 
the potential also behaves similar to the one obtained at 
the chemical potential, namely the inverse square root 
form, which indicates that the ISs are considerably nar- 
row. Also from the density profile calculated (not shown 
here, however, the results of a similar calculation can be 
found in Ref. [25|), we see that, due to the strong mod- 
ulation, the outermost IS is narrow and it's effect is neg- 
ligible, thus the slope remains almost insensitive to the 
change in the B field on this scale. An interesting com- 
parison concerning the sample widths reveals that the 
narrower the sample is, the stronger the slope. Hence 
in the Mach-Zehnder experiments and also considering 
the fact that the measurements are performed at an in- 
termediate magnetic field strength (B ~ 2.5 — 4.5T) and 
narrow samples {d ~ 1/xto), the assumption of a constant 
velocity independent of B is not realistic. We observe 
that, the disorder potential does not effect this general 
behavior as long as the dominating scattering processes 
come from the edges of the sample. Introducing disor- 
der obviously results in density fluctuations, which can 
be screened by the 2DES if the system is compressible 
(far from integer filling factors) and the conclusion is the 
opposite if the Landau levels are fully occupied. We con- 
sider a situation such that the magnetic field is tuned to 
an interval, where the average filling factor becomes close 
to an (even) integer. In this situation, a large IS is formed 
at the bulk (without modulation), spl it into several rib- 
bons (as observed in Fig. 2 of Ref. [22|) and the effect of 
these incompressible ribbons on the slope at the edge is 
marginal. This is seen in the left panel of Fig.[5j where we 
examine the behavior of the derivative comparing Vq = 
and Vq ^ 0. In the unmodulated case the slope drops 
linearly with increasing B, until a large IS is formed at 
the bulk (e.g. in Fig. [If B ~ 7.3T). The wide strip disap- 
pears when the magnetic field strength is strong enough 
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so that the Fermi level is pinned to the lowest Landau 
level, B ~ 7.45T. For the modulated case the derivative 
decreases also linearly, with a smaller slope. However, 
this linear region is larger compared to the unmodulated 
case, e.g. in figure [5]i up to B ~ 7.8T for Vb = 100 meV 
and B ~ 9.0 T for Vb = 200 meV. Depending strongly on 
the modulation amplitude, the rapid decay of the slope 
due to the formation of a large bulk IS, is observed in a 
relatively narrow B interval. The "linear slope regime" 
is observed for all considered sample widths; however, for 
narrow samples the B interval is larger for higher modu- 
lation amplitudes. This indicates that for high mobility 
samples, where the long-range part of the disorder poten- 
tial is well screened 2 ^, the linear regime will be observed 
in a narrow B interval. From the above discussion we 
conclude that, the electron velocity on the ISs presumes 
a linear B field dependence. At this point we find it useful 
to make a connection between our results and the Mach- 
Zehnder type samples. These samples have intermediate 
mobility and are relatively narrow. We have shown that 
the electron velocity calculated at the chemical poten- 
tial decreases like an inverse square root of the B field 
and the assumption of constant velocity is not applicable. 
Whereas, if the current is carried by the ISs, assuming a 
constant v e i in the magnetic field interval where the inter- 
ference pattern is observed is still irrelevant. Recently it 
has been shown theoretically that2S, within the screening 
picture of integer quantized Hall effect, the interference 
can be observed only in a narrow magnetic field inter- 
val within the plateau regime at high mobilities. The 
boundaries to observe interference pattern is estimated 
such that two separated ISs should be formed (similar 
to B < 8T of figure [Bf ) which are larger than the Fermi 
wave length (B > 6T). This interval coincides with the 
linear velocity regime shown above. Therefore we sup- 
port the idea^, that the phase of the electron calculated 
within the single particle picture should also be recon- 
sidered form the interaction point of view as presented 
in this work. So far we have examined the magnetic 
field dependence of the slope of the screened potential 
at the Fermi level and within the ISs. We have found 
that, depending on the B strength, the electron veloc- 
ity exhibits different behaviors depending on where the 
slope is calculated. If it is assumed that the current flows 
from the Landauer-Buttiker type edge states the velocity 
takes the form B^ 1 / 2 . Whereas, if the current is carried 
by the ISs, the velocity drops linearly in the case of two 
separate ISs and is highly non-linear in the presence of 
a large IS in the bulk. Next, we discuss the extend of 
the ISs depending on the magnetic field and steepness 
of the confinement potential considering different sample 
widths and temperatures. 



V. THE FORMATION OF IS 

The long-standing question of "where the current flows 
?" in the quantum Hall bar systems has been addressed 



in many different theoretical work o 10 ' 11 ^ 27 ' 28 ! 29 ' 30 . In a 
recent model using a local version of the Ohm's law, it 
was shown that the external current is confined in the ISs 
where the longitudinal resistivity vanishes, i.e. pi(x) = 0. 
This novel approach brought a quantitative explanation 
to many interesting aspects of the integer quantized Hall 
effect, among which are the high reproducibility of the 
very accurate quantized Hall plateaus, the transition be- 
tween the zero states and the description of the local cur- 
rent distribution. This model is based on the formation 
(and disappearance) of the ISs, and now we concentrate 
on their widths taking into account different edge profiles 
and sample widths. 

In Fig. O we plot the widths of the ISs (W2, for lo- 
cal filling factor 2) against the magnetic field strength 
considering different sample properties. For the constant 
donor distribution (c = 1) we see that the sample width 
has no influence on the B dependence of the width of 
the strips. Whereas W% increases by increasing the sam- 
ple width, as expected. Note that, since the variation of 
the self-consistent potential at the IS is hio c , the slope 
is calculated simply by dividing this variation by Wi- 
As a direct consequence, the slope becomes small when 
W2 becomes large. At the first glance, for intermediate 
steepness (c = 0.8) the functional form of the inhomo- 
geneous donor distribution, i.e. pi(x) or P2(x), has no 
influence on IV2. However, the linear velocity regime is 
much more extended for the etched samples than of the 
doped ones for relatively large sample widths (d > 3/im). 
This feature is more pronounced for c = 0.6, i.e. for 
the steeper edge profile, and the linear (velocity) regime 
is observed in a larger magnetic field interval compared 
to other steepness parameters. The slope of the linear 
regime is smaller for the etched sample and a smoother 
transition to the non-linear regime is observed for the 
doped edge profile, whereas the functional form of the 
donor distribution seems to show no important difference 
for different sample widths. 

As a final remark on the IS widths, we would like to 
recall the findings of Giiven and Gerhardts^ where the 
high current regime was also investigated. It was shown 
that, a large imposed current leads to a broadening of the 
ISs on one side of the sample, hence a change in the slope, 
which was also supported by the experiments^. This re- 
sult shows that there is a relation between the amplitude 
of the imposed current and the average electron velocity 
inside the ISs. We believe that, the investigation of the 
out-of-the-linear response regime will improve our under- 
standing of the Mach-Zehnder type of interferometer ex- 
periment. Our preliminary results show that, the widths 
of the ISs increase linearly by increasing the amplitude 
of the applied current. 



VI. SUMMARY 

In summary, we have calculated the slope of the self- 
consistent potential, within the Thomas-Fermi-Poisson 



8 



theory of screening. We considered two different pictures 
of edge states, namely the single particle and the incom- 
pressible states, to obtain electron velocities in the pres- 
ence of a strong perpendicular magnetic field. We have 
systematically investigated the effect of the sample prop- 
erties such as the sample width, edge profile and disorder 
potential on the electron velocities. 

We first obtained a functional form of the v e \ and the 
electric field depending on the magnetic field strength, 
without taking into account the formation of the ISs, 
and considering only the Landau quantization. It is 
shown that the interpretation of the SPV experiments 1 ^ 
strongly contradicts with our results and also with other 
experiments 19 . We found that the slope of the self- 
consistent potential changes as ~ £U 1,/2 , whereas the 
electric field at the edge behaves as F oc B 1 ^ 2 . We also 
concluded that assuming a constant v e i may lead to dis- 
crepancies in analyzing the results of Mach-Zehnder in- 
terferomete r 1 ' 15 type experiments. 

Secondly, by evaluating the full self-consistent poten- 
tial, we were able to obtain the electron velocities at the 
Fermi level and within the ISs. We found that, the full 
self-consistent results coincide with our semi-consistent 
findings pointing to the inverse square root dependence 
of v c \ within the single particle picture. The slope of the 



fully screened potential calculated at the ISs, however, 
exhibits two different regimes of magnetic field. These 
two regimes are identified by the dependence of the elec- 
tron velocity on the magnetic field which is linear in one 
regime and non-linear in the other. 

Our results indicate that, in narrow Hall bar geome- 
tries with intermediate mobilities, the edge profile be- 
comes very important in determining the electron veloc- 
ity for both the Landauer-Buttiker or the IS type edge 
states. It appears to us that re-examining the results of 
Mach-Zehnder interferomete r 14 ' 15 experiments from self- 
consistent point of view will thus be helpful to understand 
the underlying physics of the obtained interference pat- 
terns. 
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